Non-linear conductivity and quantum interference in disordered metals 
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We report on a novel non-linear electric field effect in the conductivity of disordered conductors. 
We find that an electric field gives rise to dephasing in the particle-hole channel, which depresses 
the interference effects due to disorder and interaction and leads to a non-linear conductivity. This 
non-linear effect introduces a field dependent temperature scale Te and provides a microscopic 
mechanism for electric field scaling at the metal-insulator transition. We also study the magnetic 
field dependence of the non-linear conductivity and suggest possible ways to experimentally verify 
our predictions. These effects offer a new probe to test the role of quantum interference at the 
metal-insulator transition in disordered conductors. 

PACS numbers: 72.10-d, 72.15.Rn 



Disordered conductors have been the subject of theo- 
retical and experimental study for almost twenty years 
Recently there has been a strong resurgence of 
interest in the field due to the unexpected discovery of 
a metal-insulator transition in two-dimensional systems 
H . Various suggestions have been made concerning the 
origin of the temperature dependence of the resistivity in 
the metallic phase and the nature of the metal-insulator 
transition @. One main question is whether the transi- 
tion is of a classical origin or if it is a real quantum phase 
transition. In the first case, if a standard Landau quasi- 
particle picture applies the observed resistivity could be 
attributed to a temperature dependent scattering time in 
the context of the semi-classical Boltzmann-Landau ki- 
netic equation |j| . In the second case, it has been pointed 
out H that the occurrance of a metallic phase and a 
metal-insulator transition in two dimensional systems is 
indeed possible within the standard theory of disordered- 
interacting electrons . 

To discriminate between these possibilities one needs 
specific probes for quantum interference effects. Weak lo- 
calization (WL) is probed effectively by the application 
of a magnetic field. This magnetic field also affects the 
quantum interference from the combined contribution of 
disorder and electron-electron interaction (EEI) || but 
only in the particle-hole triplet channel due to the Zee- 
man coupling. In this paper we propose a new probe for 
the EEI contribution (both the singlet and triplet) based 
on the non-linear conductivity in the presence of a static 
electric field. We recall that WL is not affected by such 
afield |g. 

To be more specific we have found that a static (or low 
frequency) electric field: (1) acts as a source of dephasing 
in the particle-hole channel and introduces a character- 
istic temperature Te = (D qp e 2 E 2 )^ below which inter- 
ference effects are suppressed. Here D qp is the quasi- 
particle diffusion coefficient defined below. (2) For tem- 



peratures above Te non-linear effects in the conductivity 
appear as T E /T 3 corrections. (3) For large electric field 
the scale of the magnetic field in the magnetoresistance 
is set by Te- Clearly one expects that besides quantum 
interference, heating effects will also be important in the 
non-linear conductivity and we shall suggest how to dis- 
tinguish the two effects. Our results, besides providing a 
new probe for EEI corrections, offer a microscopic mech- 
anism for the electric field scaling which is observed in 
two-dimensional systems |T0| . 

Before giving details of the mathematical derivation a 
qualitative understanding of the effect may be obtained 
by simple physical arguments along the lines of Ref . jll] . 
In a generalised Hartree-Fock picture one electron is scat- 
tered by the potential created by all the other electrons. 
Due to disorder, the electron density is not uniform and 
hence this potential is random. From a semi-classical 
point-of-view, a local, single-particle quantity, like cur- 
rent, only involves closed paths. Futhcrmore, the cor- 
rections are dominated by all the other electrons retrac- 
ing backwards-in-time (as holes) almost the same closed 
paths. According to Ref. [jll| only trajectories which are 
traversed in a time rj < 1/T contribute to quantum cor- 
rections. Although the two electrons go around the same 
closed path they have different starting positions. The 
first electron starts at the observation point Xi at time 
zero, while the second electron will only start to retrace 
the path at the point of interaction x 2 at time t%. This 
means that the second electron will traverse part of the 
closed path at a different time. In the presence of a vec- 
tor potential the accumulated phase difference is then 
01 - 02 = e J^dfixff) • A(x x (*'),*') -e\l df x 2 (f) • 
A(x2(i')> t'). If the vector potential is time independent 
(eg. a magnetic field) these phases completely cancel. 
However, if the vector potential is time dependent (as 
for a static electric field) the time delay leads to a finite 
phase difference 4>i — 02 = e(x 2 — xi) • E77, which suggests 
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that the EEI correction should be sensitive to a static 
electric field, in contrast to WL. Such a phase-sensitivity 
leads to non-linear conductivity. It is possible to esti- 
mate the typical electric field scale where dephasing and 
non-linear effects in a weakly disordered metal become 
important. The typical time scale is the inverse temper- 
ature and the typical length scale is the thermal length 
Lt = (Dqp/T) 1 / 2 . The non-linear effects become impor- 
tant when the phase difference induced by the electric 
field is of order one, which leads to the condition that 
the voltage drop over a thermal length becomes of the 
order of the temperature, i.e., when eELx ~ T. This 
condition defines the temperature scale Tg given above. 

We now present a quantitative theory of our results. 
We start with the expression for the EEI quantum correc- 
tion to the current due to the interplay between disorder 
and interaction. Within the real-time Keldysh formalism 
we obtain: 
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A pictorial representation of this equation is shown in 
Fig.|l|. The details of its derivation may be found in jl3| . 
The sum ^] JM is over one singlet (J = 0, M = 0) and 
three triplet channels with J = l,M = 0, ±1. In eq.(|l]) 
r is the elastic scattering time which is the short-time 
cut-off in the problem. Vj t M and D r j M (t,t') are the in- 
teraction and the diffusion propagator in the spin channel 
(J, M). Here the time arguments t, t' refer to the incom- 
ing and outgoing centre-of-mass time of the particle-hole 
pair and 77 to the relative time which is constant during 
the propagation. Notice that both Vj : m and -Dj.m are 
retarded functions. The factor containing the hyperbolic 
sine comes from the Fourier transform of Fermi functions 
and limits us to trajectories with traverse time 77 < 1/T. 
The interaction is found by summing ladder diagrams 
and is given by 



Vj,Af(q,w) = 7,7 — 
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where 77 is the static amplitude for scattering between 
particles and holes. The quasi-particle diffusion constant 
can be expressed in terms of the particle diffusion con- 
stant -Dqp = D/Z. The parameter Z, which arises in the 
context of the Fermi liquid theory of disordered systems 
as the energy renormahzation, plays the role of mass 
renormahzation, m*/m, in the effective Fermi liquid the- 
ory of disordered systems 114]. The interaction ampli- 
tude in the spin singlet channel is given by 77=0 = 1/2 
for long range Coulomb forces. The triplet amplitude, 



for which we adopt in the following the standard nota- 
tion 7,7=1 = —72/2, is related to the Landau parameter 
F° via 72 = -A° a = -F Q °/(1 + F°), The diffusion prop- 
agator Dj.m is given by the solution of the differential 
equation 
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where A, ; (r,/j) = A(r,t + rj/2) - A(r,< - 77/2). The 
term iMfl s is due to the Zeeman coupling, where fl s = 
(1 + 72)^5 with £l s = g/isH. 

For a better understanding of eq.([j]) we make contact 
with the physical arguments given in the introduction. 
By considering Fig.[j] one observes that for a piece of the 
path the particle and hole are delayed by a time 77, the 
traverse time of the closed path. This corresponds to 
the second of the diffusons in eq. (Q) . In the other piece 
there is no delay between the particle and hole. This 
corresponds to the first diffuson in eq.(|]). 

We now evaluate the current explicitely. According to 
eq.(H) the interaction Vj,jvf and the first of the two dif- 
fusons in (|l|) do not depend on the vector potential. An 
electric field, however, affects the second diffuson in (|l|) 
due to the non-zero time delay 77 between the particle and 
hole. For a static field the vector potential is A(t) = — Et 
and the solution of (pi) is D v , M (t — | , t\ — 
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The 
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equation for the current after integrating over the mo- 
mentum then becomes 
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where we have introduced T% = D qp (eE) 2 and d is the 
dimension. It is clear from this equation that the electric 
field provides a dephasing time ~ , since in the low 
temperature limit T <C Te the exponential now cuts off 
all times larger than T^ 1 . 

We first consider the current in the weak electric field 
regime and derive the leading non- linear terms. In the 
absence of magnetic field we find 
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where the functions /^' 3 are shown in the table. For the 
sake of completeness we have also included the term lin- 
ear in the electric field which reproduces the well-known 
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interaction correction to the conductivity. Notice that 
the functions f^' 3 ^) are the sum of the singlet and the 
triplet contributions. The remaining integrals are of a 
standard form. We now have the results for 5a = <5|j|/|E| 



5a 2 = 
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where Lt is the thermal length defined previously and 
we have left out temperature independent terms. We 
recall that, in the case of spin-singlet interactions only 
(72 = 0) the conductivity decreases with decreasing tem- 
perature (i.e. /J(0) > 0), whereas for sufficiently large 
triplet amplitude 72 the latter dominates and leads to an 
increase of conductivity with decreasing temperature (i.e. 
/K72) < 0). The non-linear coefficient /J (72), however, 
is generically positive and only changes sign for large 72 
in d = 3. 

We study the cross-over behaviour from small to large 
electric fields by numerically integrating eq.(^|). The con- 
ductivity as a function of the electric field is plotted 
in figj| for two values of 72 for d = 2. At zero field, 
for 72 = (72 = 5) the correction is localising (anti- 
localising) with 5a 2 < (5a 2 > 0). The quadratic in- 
crease at small fields has a positive curvature irrespective 
of the value of 72 . At large field, the temperature scale 
disappears and the correction 5a has the same form as 
the linear conductivity with Te replacing the tempera- 
ture. 

Non-linear effects also appear in the magnetoconduc- 
tance which originate from the magnetic field depression 
of the M = ±1 triplet contributions to the current. In 
particular we find for small Te and small Zeeman energy 
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where Aa = cr(O s ) — cr(0) and the functions g^ 3 are 
also shown in the table |lq] . To illustrate the behaviour 
at large Q s we again resort to numerical integration of 
eq.(Q). In figj^ we show the magnetic field dependence 
of the current for two choices of Tg for 72 = 2.5. Notice 
that for such a value of 72 the zero magnetic field conduc- 
tivity interference correction has an anti-localizing char- 
acter. This explains the suppression of conductivity with 
increasing Tg. For small Te we obtained the standard 
behaviour of a initial quadratic decrease on the scale of 



the temperature followed by a logarithmic suppression of 
the corrections. For large Tg however, the temperature 
disappears as an energy scale and, although the curve ap- 
pears similar, the scale of the magnetic field is now set by 
Te- Expanding eq.(Q) to leading order in the magnetic 
field for T E -»T one obtains Act cx ft 2 /T 2 . 

The effects described in this paper may be detected by 
measuring the current-voltage characteristics. In such a 
measurement however the electron temperature changes 
with the applied voltage and one has to discriminate 
heating from non-heating non- linear effects. A direct way 
to isolate the non-linear contribution due to the dephas- 
ing effect of E would be to measure the electron tempera- 
ture T e i for a given E (for instance by noise measurements 
as in @). Then <r(T e/ ,0) - a(T,E) yields the effect of 
the electric field on the EEI contribution and provides 
a direct probe of the relevance of quantum interference 
in the p-h channels. Alternatively, at low temperature, 
where T e \T e \- p h 1 (r e i_ p h is the heat electron-phonon 
relaxation time) non-heating effects could be detected by 
exploiting the different time scales r e i_ p h and T^ 1 which 
control the frequency dependence of heating and non- 
heating effects respectively. In a time-dependent elec- 
tric field, E(t) — Ecos(ftt), the electron temperature 
becomes time dependent. For frequencies £1 > 1/T c i_ p h 
however the temperature cannot follow the electric field, 
i.e. heating becomes time independent. Non-linearities 
due to quantum interference on the other hand follow 
the electric field instantly as long as the frequency re- 
mains smaller than the temperature. Thus measuring 
non-linear response in the presence of a microwave with 
a frequency of the order CI > 1/T c i_ p h offers a possibility 
to detect the predicted effects p7| . 

Possible non-heating effects have already been ob- 
served in different materials like 2D Si-MOSFETs @,||, 
GaAs @, SiGe §o| and in gold films §l],||. In partic- 
ular in Rcf. jl(| a remarkable electric field scaling was 
observed near the metal-insulator transition. Our theory 
of the non-linear effects provides an explicit mechanism 
for the electric field scaling via the temperature scale Te'- 
on the basis of general scaling arguments close to a quan- 
tum critical point, the temperature scales as T ~ £;~ z 
where £ is the correlation length and z is the dynamical 
critical exponent |p3f . In a diffusive system temperature 
and length scales are related via the diffusion constant 
with T ~ -D qp (£)/£ 2 implying a scaling of D/Z = D qp 
near the critical point as D qp ~ £ 2 ~ z . From our rela- 
tion Tj = L> qp e 2 £; 2 one then obtains E - ^ ( - 1+z \ In 
the experiments z is near one, which corresponds to a 
growing quasi-particlc diffusion constant and a vanishing 
quasi-particle density of states near the transition. This 
small value of z < 2 implies that the electronic specific 
heat would vanish as c„ ~ T£ z ~ 2 ~ T 2 / z . From these 
considerations one expects large non-linear effects near 
the critical point. Such large effects have been observed 
in Ref. |HJ for a GaAs metallic sample near the metal- 



3 



insulator transition. In this experiment the differential 
conductivity was found to decrease with increasing volt- 
age. In particular at low temperature the conductivity 
shows a non-monotonic behaviour of the type shown in 
fig. 2 for large 72. By comparing with experiment it ap- 
pears that the scale over which the effect is observed is 
larger than what we would predict based upon a naive 
estimate of the diffusion constant from the conductiv- 
ity a — 2e 2 DN and the free particle density of states 
Nq (i.e. assuming D qp ~ D). However, by allowing for 
quasi-particle diffusion constant renormalization as im- 
plied by scaling one can obtain larger effects. 

Finally, we point out that the effects discussed in this 
paper could have an enhanced relevance in the presence 
of strong local electric fields such as in percolative metal- 
lic systems. 
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FIG. 1. Pictorial representation of the current formula. 
Four dashed lines represent a diffuson. 




FIG. 2. The electric field dependence of the interaction cor- 
rection to the conductivity in two dimensions in units of e 2 /h 
for 72 = (dashed line) and 72 = 5 (solid line). Note the 
different scales used for the two values of 72. 
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FIG. 3. The magnetic field dependence of the conductiv- 
ity in two dimensions in units of e 2 /h for two different val- 
ues of Te and for 72 = 2.5. The inset shows the difference 
AT E a2(O s ) = <5<J2(5T, Q s ) — 8a2(0.1T, fi s )of the two curves. 
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TABLE I. Table of coefficients which appear in the ex- 
pressions for the current. For small 72 these reduce to 
/K72) = 2/d - 372/2, f! (72) = 4/(d(2 + d)) - 72/4, 
3d(72) = 72/2 and ^(72) = -72/12. 
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